Burst Mechanisms in Hydrodynamics 

E. Knobloch & J. Moehlis 

Department of Physics, University of California, 

Berkeley, CA 94720, USA 

E-Mail: knobloch@physics. berkeley. edu 

Abstract 

Different mechanisms believed to be responsible for the generation of bursts 
in hydrodynamical systems are reviewed and a new mechanism capable of 
generating regular or irregular bursts of large dynamic range near threshold 
is described. The new mechanism is present in the interaction between os- 
cillatory modes of odd and even parity in systems of large but finite aspect 
ratio, and provides an explanation for the bursting behavior observed in 
binary fluid convection. Additional applications of the new mechanism are 
proposed. 

1 Introduction 

Bursts of activity, be they regular or irregular, are a common oc- 
currence in physical and biological systems. In recent years several 
models of bursting behavior in hydrodynamical systems have been 
described using ideas from dynamical systems theory. In this arti- 
cle we review these and then describe a new mechanism (Moehlis &: 
Knobloch ^ , Knobloch & Moehlis @ ) which provides an explana- 
tion for the bursting behavior observed in experiments on convection 
in ^He/^He mixtures by Sullivan & Ahlers [^. This mechanism ope- 
rates naturally in systems with broken D4 symmetry undergoing a 



Hopf bifurcation from a trivial state. This symmetry, the symmetry 
group of a square, may be present because of the geometry of the 
system under consideration (for example, the shape of the container) 
but also appears in large aspect ratio systems with reflection sym- 
metry (Landsberg & Knobloch [^). In either case bursting arises as 
a result of the nonlinear interaction between two nearly degenerate 
modes with different symmetries. 

This article is an expanded version of Knobloch Sz Moehlis . 

2 Mechanisms producing bursting 

As detailed further below, bursts come in many different forms, dis- 
tinguished by their dynamic range (i.e., the range of amplitudes dur- 
ing each burst), duration, and recurrence properties. Particularly 
important for the purposes of the present article is the question of 
whether the observed bursts occur close to the threshold of a primary 
instability or whether they are found far from threshold. In the for- 
mer case a dynamical systems approach is likely to be successful: in 
this regime the spatial structure usually resembles the eigenfunctions 
of the linear problem and it is likely that only a small number of 
degrees of freedom participate in the burst. Such bursts take place 
fundamentally in the time domain with their spatial manifestation 
of secondary importance, in contrast to pulses which are structures 
localized in both time and space; the latter are not considered here. 
Since the equations governing the evolution of primary instabilities 
are often highly symmetric (see Crawford & Knobloch Q) global bi- 
furcations are likely to occur and these serve as potential candidates 
for bursting mechanisms. In contrast, bursts found far from threshold 
usually involve many degrees of freedom but even here some progress 
is sometimes possible. 

2.1 Bursts in the wall region of a turbulent boundary 
layer 

The presence of coherent structures in a turbulent boundary layer is 
well established (see, e.g., Robinson Q and the collection of articles 



edited by Panton Q). The space-time evolution of these structures 
is often characterized by intermittent bursting events involving low 
speed streamwise "streaks" of fluid. Specifically, let xi,X2, and X3 be 
the streamwise, wall normal, and spanwise directions with associated 
velocity components U + ui,U2, and ^3, respectively; here U{x2) is 
the mean flow. In a "burst" the streak breaks up and low speed fluid 
moves upward away from the wall (ui < 0,^2 > 0); this is followed 
by a "sweep" in which fast fluid moves downward towards the wall 
{ui > 0,U2 < 0). After the burst/sweep cycle the streak reforms, 
often with a lateral spanwise shift. 

A low-dimensional model of the burst/sweep cycle was developed 
by Aubry et al. [0]; further details and later references may be found 
in Holmes et al. [11, 12|. To construct such a model the authors used 
a Karhunen-Loeve decomposition of the data to identify an energet- 
ically dominant empirical set of eigenfunctions, hereafter "modes". 
The original study of Aubry et al. [10| used experimental data for 
pipe flow with Re ~ 6750, while later studies used data for chan- 
nel flow from direct numerical simulation with Re ~ 3000 — 4000 
and large eddy simulations with Re ^ 13800 (Holmes et al. |12]). 
The model was constructed by projecting the Navier-Stokes equation 
onto this basis and consists of a set of coupled ODEs for the ampli- 
tudes of these modes. The fixed points of these equations are to be 
associated with the presence of coherent structures. There are two 
types, related by half-wavelength spanwise translation. 

Numerical integration of the model reveals that these fixed points 
are typically unstable and that they are connected by a heteroclinic 
cycle. In such a cycle the trajectory alternately visits the vicinities 
of the two unstable fixed points. In the model of Aubry et al. |[l^ 
this heteroclinic cycle is found to be structurally stable, i.e., it per- 
sists over a range of parameter values. This is a consequence of the 
0(2) symmetry of the equations inherited from periodic boundary 
conditions in the spanwise direction. Moreover, for the parameter 
values of interest this cycle is attracting, i.e., it attracts all nearby 
trajectories. Since the transition from one fixed point to the other 
corresponds to a spanwise translation by half a wavelength, the re- 



current excursions along such a heteroclinic cycle can be identified 
with the burst/sweep cycle described above. However, since this cy- 
cle is attracting, the time between successive bursts will increase as 
time progresses. This is not observed and Aubry et al. [10| appeal 
to the presence of a random pressure term modeling the effect of the 
outer fluid layer to kick the trajectory from the heteroclinic cycle. In 
the language of Busse |13| such a pressure term results in a statisti- 
cal limit cycle, with the bursting events occurring randomly in time 
but with a well-defined mean rate. The resulting temporal distribu- 
tion of the burst events is characterized by a strong exponential tail, 
matching experimental observations. 

Attracting structurally stable heteroclinic cycles occur in a num- 
ber of problems of this type, i.e., mode interaction problems with 
0(2) symmetry (Armbruster et al. [14|, Proctor & Jones |15], Mel- 
bourne et al. H^, Steindl & Troger |17], Krupa [18|, Hirschberg & 
Knobloch H]). 



2.2 Bursts in shear flows undergoing subcritical tran- 
sition to turbulence 

Experimental studies of plane Couette flow and Poiseuille (pipe) flow 
have shown that at high enough values of the Reynolds number Re the 
basic laminar flow becomes turbulent. However, these laminar flows 
are linearly stable at all Re (see, e.g., Drazin & Reid (2^]). Con- 
sequently, the transition to turbulence in these systems must arise 
from finite (i.e., not infinitesimal) perturbations to the basic flow; 
such transitions are called subcritical since the turbulent state exists 
for values of Re for which the laminar state is stable. Much recent 
work (reviewed in Baggett & Trefethen ||2l|) has emphasized the im- 
portance of the fact that the linear operators L governing the evolu- 
tion of perturbations of these flows are non-normal, i.e., L^L ^ LL\ 
where is the adjoint of L. Linear systems with a non- normal L 
can exhibit transient growth even though the laminar state is linearly 
stable; if the growth is large enough, nonlinearities in the system may 
then trigger a transition to turbulence. Analysis of low-dimensional 
models supported by numerical simulations suggests that the mini- 



mum pertm'bation amplitude e that results in turbulence scales as 
e = 0(i?e°) for some q < -1 (Baggett & Trefethen Thus, e 

decreases rapidly with increasing Re, and the experimentally deter- 
mined Re for transition should be that value at which e is roughly 
equal to perturbations due to noise or imperfections in the system. 

The turbulence excited by finite amplitude perturbations in shear 
flows often takes the form of turbulent spots which can move, grow, 
split, and merge (see, e.g., Daviaud et al. ||2^). Turbulent spots can 
also burst intermittently. For example, Bottin et al. |23] observed 
intermittent turbulent bursts in plane Couette flow for Re ~ 325; 
these bursts were triggered by a spanwise wire through the center of 
the channel, and could be localized in space by introducing a bead 
into the central plane. We focus here on this burst-like behavior, also 
found in low-dimensional models of the subcritical transition to tur- 
bulence. The important issue in such models is the attractor to which 
the flow settles after an appropriate finite perturbation to the basic 
laminar flow. This will depend crucially on nonlinear terms in the 
equations, a point emphasized, for example, by Waleffe ^4, 25, 27] 
and Dauchot & Manneville |2^. If this attractor is a fixed point, the 
system settles into a steady state in which bursts do not occur. On 
the other hand, if this attractor is a limit cycle it may be appro- 
priate to interpret the resulting behavior as burst-like. This is the 
case in the model studied in Waleffe [24| in which the amplitudes of 
streamwise streaks, streamwise rolls, and the streak instability can 
periodically undergo short-lived explosive growths at the expense of 



the mean shear. Waleffe |24, pq| refers to this as a self-sustaining 



process. Since the laminar state is linearly stable, such limit cycles 
cannot bifurcate off the laminar state. Instead, as pointed out in 
Waleffe |2^, |2^, they may be born in a Hopf bifurcation (from a 
fixed point other than that corresponding to the laminar state) or 
a homoclinic bifurcation in which a trajectory homoclinic to a fixed 
point forms at some value of Re such that for smaller (larger) values 
of Re a limit cycle exists (does not exist), or vice versa. Near such 
a homoclinic bifurcation, bursts are expected because the trajectory 
spends a long time near the fixed point, bursting away and return- 



ing in a periodic or chaotic fashion, depending on the eigenvalues of 
the fixed point. Other generic codimension one bifurcations leading 
to the appearance or disappearance of limit cycles are saddle-node 
bifurcations of limit cycles and saddle-nodes in which a pair of fixed 
points appears on the limit cycle (Guckenheimer et al. ||2^). Limit 
cycles could also come in from or go off to infinity. In other models 
chaotic attractors have been found (see, e.g., Gebhardt & Grossman 
3C] and Baggett et al. ||3l|) and these may also give rise to burst-like 



behavior. 

The relation of the bursts described in this section (which occur 
for moderate values of Re, e.g.. Re ~ 325) to those which occur in 
the turbulent boundary layer (which occur for high values of Re) is 
not completely clear. However, the similarities in the phenomenology 



have led Waleffe [27] to propose that the self-sustaining process will 



continue to have importance in the near-wall region of high Re flows. 

2.3 Heteroclinic connections to infinity 

Another mechanism involving heteroclinic connections, distinct from 



that discussed in section 2A, has been investigated by Newell et 
al. 1 32, 33 1 as a possible model for spatio-temporal intermittency 
in turbulent flow. The authors suggest that such systems may be 
viewed as nearly Hamiltonian except during periods of localized in- 
tense dissipation. A related "punctuated Hamiltonian" approach to 
the evolution of two-dimensional turbulence has met with consider- 



able success (Carnevale et al. |34] and Weiss & Mc Williams ]35[|). 



For their description Newell et al. divide the instantaneous states of 
the flow into two categories, a turbulent soup (TS) characterized by 
weak coherence, and a singular (S) state characterized by strong co- 
herence, and suppose that the TS and S states are generalized saddles 
in an appropriate phase space. Furthermore, they suppose that in the 
Hamiltonian limit the unstable manifold of TS (S) intersects transver- 
sally the stable manifold of S (TS). If the constant energy surfaces are 
noncompact (i.e., unbounded), the evolution of the Hamiltonian sys- 
tem may take the system into regions of phase space with very high 
("infinite") velocities and small scales. These regions are identified 



with the S states and high dissipation. In such a scenario the strong 
dissipation events are therefore identified with excursions along hete- 
rochnic connections to infinity. Perturbations to the system (such as 
the addition of dissipative processes) may prevent the trajectory from 
actually reaching infinity, but this underlying unperturbed structure 
implies that large excursions are still possible. 

Newell et al. apply these ideas to the two-dimensional nonlin- 
ear Schrodinger equation (NLSE) with perturbations in the form of 
special driving and dissipative terms which act at large and small 
scales, respectively. Here S consists of "filament" solutions to the un- 
perturbed NLSE which become singular in finite time and represent 
coherent structures which may occur at any position in the flow field. 
When the solution is near S a large portion of the energy is in small 
scales; for the perturbed equations the dissipative term then becomes 
important so that the filament solution is approached but collapses 
before it is reached. This leads to a spatially and temporally random 
occurrence of localized burst-like events for the perturbed equation. 
The rate of attraction at S is determined by the faster than expo- 
nential rate at which the filament becomes singular, while the rate 
of repulsion at S is governed by the dissipative process and hence is 
unrelated to the rate of attraction. 

This bursting mechanism shares characteristics with that de- 
scribed in Kaplan et al. |36] in which solutions of a single complex 
Ginzburg-Landau equation with periodic boundary conditions un- 
dergo faster than exponential bursting due to a destabilizing non- 
linearity and collapse due to strong nonlinear dispersion (see also 
Bretherton & Spiegel [^). A study of a generalization of Burger's 
equation modeling nonlocality effects suggests the presence of burst- 
like events through a similar scenario (Makarenko et al. Isl 



2.4 Bursts in the Kolmogorov flow 

The Kolmogorov flow u = (A; sin /cy, 0) is an exact solution of the 
two-dimensional incompressible Navier-Stokes equation with unidi- 
rectional forcing f at wavenumber k: i = {uk^ sinky,0). With in- 
creasing Reynolds number Re = z/~^ this flow becomes unstable, 



and direct numerical simulation with 27r-periodic boundary condi- 
tions shows that for moderately high Reynolds numbers and k > 1 
the resulting flow is characterized by intermittent bursting (She |39|] , 



Nicolaenko & She [|g, |2|, Armbruster et al. @ ||, 
A burst occurs when the system evolves from a coherent vortex-like 
modulated traveling wave (MW) to a spatially disordered state fol- 
lowing transfer of energy from large to small scales. The system then 
relaxes to the vicinity of another symmetry-related MW state, and 
the process continues with bursts occurring irregularly but with a 
well-defined mean period. 

The details of what actually happens appear to depend on the 
value of k because k determines the symmetry of the nonlinear equa- 
tion describing the evolution of the perturbation streamfunction (p 
about the Kolmogorov flow. Although there is no compelling reason 
for it, all simulations of this equation have been performed with 27r- 
periodic boundary conditions in both directions. With these boun- 
dary conditions this equation has a symmetry which is the semi- 
direct product of the dihedral group D2k (generated by the actions 
{x,y,(j))^ {-x,-y,(t)) and {x,y,cl)) {-x,y + tt /k, -(p)) and the 
group S0(2) (representing the symmetry under translations x — > 
x+const). In the simplest case, k = 1, this symmetry group is iso- 
morphic to the direct product of the group 0(2) of rotations and 
reflections of a circle and the group Z2 representing reflections in the 
y-direction. Unfortunately, when k = 1 the Kolmogorov flow with 
27r-periodic boundary conditions is not unstable for any value of Re 
(Meshalkin & Sinai , Green , Marchioro ) , and one is forced 



to consider k > 1. Armbruster et al. |46| analyzed carefully the k = 2 
case and showed that while a heteroclinic cycle between the MW 
states does form, it is not structurally stable; the bursts are therefore 
not produced by a mechanism of the type described in section 2.1. It 
is possible, however, that the onset of bursting is associated with a 
symmetry-increasing bifurcation at Re = Reg (see, e.g., Rucklidge & 
Matthews [Q). This would explain why the system stays in a single 
MW state for Re just below Reg but visits the vicinity of different 
but symmetry-related MW states for Re just above Res- However, 



despite much work a detailed understanding of the bursts in this sys- 
tem remains elusive. 

An alternative approach is to consider periodic domains with 
different periodicities in the two directions. In particular, if we con- 
sider the domain {— vr < x < vr,— vr/A; < y < n/k} with k > 1 the 
symmetry group remains 0(2) XZ2 but sufficiently long perturbations 
now grow. The unstable modes are either even or odd under the re- 
flection {x,y) {—x,—y) with respect to a suitable origin. Mode 
interaction between these two steady modes can result in a sequence 
of transitions summarized in figure I (Hirschberg & Knobloch 
the Kolmogorov flow loses stability to an even mode (Z), followed by 
a steady state bifurcation to a mixed parity state (MM^/2)- Since 
each of these states is defined to within a translation in x modulo 2tt, 
we say that it forms a circle of states. The state then loses 

stability in a further steady state bifurcation to a traveling wave 
(TW) which in turn loses stability at a Hopf bifurcation to a MW. 
The MW two-torus terminates in a collision with the two circles of 
pure parity states, forming an attracting structurally stable hetero- 
clinic cycle connecting them and their quarter-wavelength translates 
(Hirschberg &: Knobloch In this regime the behavior would re- 

semble that found in the numerical simulations, with higher modes 
kicking the system away from this cycle. Indeed this sequence of 
transitions echoes the results obtained by She |3y] and Nicolaenko & 



She Eg, 41, H|, H| for k = 8. While it is likely that the A: = 1 sce- 



nario is relevant to these calculations because of the tendency towards 
an inverse cascade in these two-dimensional systems, we emphasize 
that, as mentioned above, a careful analysis of the k = 2 case by 
Armbruster et al. [46| shows that while a heteroclinic cycle of the 
required type does indeed form, it is not structurally stable. The 
case fc = 4 has also been studied (Piatt et al. |51|) and a similar 
sequence of transitions found. Undoubtedly simulations on the do- 
main {—IT < X < vr, —Tv/k < y < ir/k} would shed new light on the 
problem, cf. Posch & Hoover fS^ . 

The group 0(2) XZ2 also arises in convection in rotating straight 
channels (Knobloch p3|) and in natural convection in a vertical slot 




Figure 1: Sample bifurcation diagram for the interaction of an odd 
and even mode for the Kolmogorov flow obtained by traversing the 
parameter space along the line 7. A structurally stable hetero- 
clinic cycle exists for a range of parameters. After Hirschberg &: 
Knobloch [19]. Other cuts through parameter space give different 
bifurcation diagrams. 

(Xin et al. []54|). In both of these cases the linear eigenfunctions are 
either even or odd with respect to a rotation by vr and the systems ex- 
hibit similar transitions. Related behavior has recently been observed 
in three-dimensional magnetoconvection with periodic boundary con- 
ditions on a square lattice. Unlike the Kolmogorov flow the evolution 
equation describing a steady state secondary instability of a square 
pattern is isotropic in the horizontal. Numerical simulations (Ruck- 
lidge et al. ||55|) show intermittent breakdown of a square pattern 
followed by its restoration modulo translation. 

2.5 Bursts in the Taylor-Couette system 

The Taylor-Couette system consists of concentric cylinders enclos- 
ing a fluid-filled annulus (see, e.g., Tagg fs^). The cylinders can be 
rotated independently. In the counterrotating regime the first state 
consists of spiral vortices of either odd or even parity with respect to 
midheight. Slightly above onset the flow resembles interpenetrating 



spiral (IPS) flow which may be intermittently interrupted by bm'sts 
of turbulence which fill the entire flow field (Hamill [5^ , Coughlin et 



al. [58, With periodic boundary conditions in the axial direction 
numerical simulations by Coughlin et al. ^] show that the IPS 
flow is temporally chaotic and consists of coexisting modes with differ- 
ent axial and azimuthal wavenumbers. This flow is confined primarily 
to the vicinity of the inner cylinder where the axisymmetric base flow 
is subject to an inviscid Rayleigh instability. Coughlin et al. |58, [59(1 
conclude that the onset of turbulence is correlated with a secondary 
instability of one of the coexisting modes of the IPS flow, namely the 
basic spiral vortex flow with azimuthal wavenumber m = 4. Indeed, 
when this mode is taken as the initial flow for parameters chosen such 
that the full IPS flow undergoes bursts, a secondary Hopf bifurcation 
from this state with the same azimuthal wavenumber but four times 
the axial wavelength is identified. The secondary instability grows 
in amplitude and ultimately provides a finite amplitude perturbation 
to the inviscidly stable flow near the outer cylinder, triggering a tur- 
bulent burst throughout the whole apparatus. During a burst small 
scales are generated throughout the apparatus leading to a rapid col- 
lapse of the turbulence and restoration of the IPS flow; the bursting 
process then repeats roughly periodically in time (see figure §). 

As discussed in section 3T, in a finite Taylor-Couette apparatus 
there is a natural mechanism for generating bursts close to onset. 
This mechanism operates in the regime e ~ T"^, where e measures 
the fractional distance above threshold for the primary instability 
and r is the aspect ratio of the annulus (Landsberg & Knobloch , 
Renardy |6^] ) . For larger e the influence of the boundaries no longer 
extends throughout the apparatus and is confined to boundary layers 
near the top and bottom. In this regime the dynamics in the bulk 
may be described by imposing periodic boundary conditions with pe- 
riod that is a multiple of the wavelength of the primary instability. 
The success of the simulation of the observed turbulent bursts using 
such periodic boundary conditions (Coughlin et al. []58| , |59| ) suggests 
that these bursts occur too far above threshold to be explained by 
the mechanism described in section 3.1 below. This suggestion is sup- 




Figure 2: Azimuthal component of vorticity in the laminar inter- 
penetrating spiral state (top), just before a burst (middle), and in 
the turbulent state (bottom), labeled using polar coordinates. After 
Coughlin et al. [59]. Courtesy K. Coughlin. 

ported by figure ^ which compares the location of the regime where 
this mechanism may be expected to operate with e;,, the experimen- 
tal value of e for the onset of bursts, as a function of T. The latter 
is obtained using the approximation ^ {Rt — Rips) / Rips-, where 
Rips and Rb are the inner cylinder Reynolds numbers for the onset 
of IPS flow and bursts, and 



Rips = (3837 ± 374)?-^ + 680 ± 14 
Rb = (3680 ± 374)r-i + 701 ± 14 
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Figure 3: The solid line shows the inner cylinder reduced Reynolds 
number ej, for the onset of bursts obtained from data in Hamill [57]. 
The dashed line shows the estimate e = near and below which the 
asymptotic bursting mechanism described in section |3.1| is expected 
to apply. The bursts found at T = 17.9 and F = 26 (Hamill [57], 
Coughlin et al. [59]) apparently fall outside the range of validity of 
the theory of section |3.l| . 

(Hamill [^^ ) . This is a good approximation for strongly counterrotat- 
ing cylinders because the IPS flow sets in for inner cylinder Reynolds 
number only about 1% above the primary onset to spiral vortices 
(Hamill |5^). The figure suggests that the observed bursts may fall 
within the range of validity of this theory for T only slightly smaller 
than those used in the experiments (F = 17.9, 26); for such bursts 
the presence of endwalls should become significant. 

2.6 Intermittency 

The term intermittency refers to occasional, non-periodic switching 
between different types of behavior; such transitions may be viewed 
as bursts. In this section we describe different types of intermittency, 
dividing them into three classes. 

2.6.1 Low-dimensional intermittency 

It is often the case that fluid dynamical systems may be modelled 
by low-dimensional Poincare maps. Fixed points of such maps corre- 
spond to limit cycles of an appropriate dynamical system; when these 




lose stability intermittency may result. Without loss of generality we 
assume that for A < Ap the system is in a "laminar" state correspond- 
ing to a stable fixed point of the map, and the onset of intermittency 
occurs at A = Ap. For A > Ap the state of the system will resemble 
the stable laminar state which existed for A < Ap for some time but 
will intermittently undergo "bursts" away from this state, followed 
by "reinjections" to the vicinity of the laminar state. The three types 
of intermittency identified by Pomeau & Manneville correspond 
to three different ways in which the laminar state ceases to exist or 
loses stability at A = Ap: Type I intermittency results when a sta- 
ble and unstable limit cycle annihilate in a saddle-node bifurcation, 
Type II results when the limit cycle loses stability in a subcritical 
Hopf bifurcation, while in Type III the limit cycle loses stability in 
an inverse period-doubling bifurcation. Type I and Type III intermit- 
tency have been observed, for example, in Rayleigh-Benard convec- 



tion (see Berge et al. ||6^ and Dubois et al. ||6^, respectively), while 
Type II intermittency has been observed in a hot wire experiment 
(Ringuet et al. ||6^). A different type of intermittency, called Type 
X, was observed by Price & Mullin [35| in a variant of the Taylor- 
Couette system; this is similar to Type I intermittency but involves 
a hysteretic transition due to the nature of the reinjection. Another 
type of intermittency, called Type V, was introduced in Bauer et al. 
6€] and He et al. |67|; this is similar in spirit to Type I intermit- 
tency but involves one-dimensional maps which are nondifferentiable 
or discontinuous. These different types of intermittency have distinct 
properties, such as the scaling behavior of the average time between 
bursts with A — Ap. 



2.6.2 Crisis-induced intermittency 

A crisis is a sudden change in a strange attractor as a parameter is 
varied (Grebogi et al. |6^, |6^, |7^). There are three types of crises, 
and without loss of generality we assume that the crisis occurs as A is 
increased through Ac. In a boundary crisis, at A = Ac the strange at- 
tractor collides with a coexisting unstable periodic orbit which lies on 
the boundary of the basin of attraction of the strange attractor. This 



leads to the destruction of the strange attractor, but chaotic tran- 
sients will still exist. Intermittency is not associated with this crisis. 
In an interior crisis the strange attractor collides with a coexisting 
unstable periodic orbit at A = Ac, but here this leads to a widening 
rather than the destruction of the strange attractor. For A slightly 
larger than Ac the trajectory stays near the region of phase space 
occupied by the strange attractor before the crisis for long times, but 
intermittently bursts into a new region. In an attractor merging crisis 
two strange attractors with basins of attraction separated by a basin 
boundary are present when A < Ac- At Ac the two strange attrac- 
tors simultaneously touch the basin boundary and "merge" to form 
a larger strange attractor for A > Ac- Such a crisis is often associated 
with a symmetry- increasing bifurcation (Chossat & Golubitsky ||7l[| ). 
For A > Ac there is a single strange attractor on which the trajec- 
tory intermittently switches between states resembling the distinct 
strange attractors which existed for A < Ac- 

2.6.3 Intermittency involving an invariant manifold 

A manifold in phase space is invariant if every initial condition on the 
manifold generates an orbit that remains on the manifold. Invariant 
manifolds often exist due to symmetries, but this is not necessary. 
There are several mechanisms for intermittency involving strange at- 
tractors on an invariant manifold. 

First, suppose that as a bifurcation parameter A is increased 
through Ap, the strange attractor on the invariant manifold loses 
stability transverse to the manifold; this is called a blowout bifur- 
cation (Ott & Sommerer |7^). Suppose that the dynamics within 
the invariant subspace do not depend on A. Such a system thus has 
a skew-product structure (Piatt et al. [^]), and A is called a nor- 
mal parameter (Ashwin et al. ^H). If the blowout bifurcation is 
"supercritical" (Ott & Sommerer |Q) then for A just above Ap a 
trajectory will spend a long time near the invariant manifold, inter- 
mittently bursting away from it, only to return due to the presence 
of a reinjection mechanism. This scenario is known as on-off inter- 
mittency, where the "off" ("on") state corresponds to the system 



being near (away from) the invariant manifold (see, e.g., Piatt et al. 
[7^ , Venkataramani et al. ||7^). Recently it has been shown that in 
appropriate circumstances a blowout bifurcation can lead to a struc- 
turally stable (possibly attracting) heteroclinic cycle between chaotic 
invariant sets (Ashwin & Rucklidge ||76[| ). 

For systems which lack a skew-product structure and thus are 
governed by non-normal parameters, in-out intermittency is possible 
(Ashwin et al. [^], Covas et al. [^]). Here the attraction and repul- 
sion to the invariant subspace are controlled by different dynamics. 
This occurs when the attractor for dynamics restricted to the in- 
variant subspace is smaller than the intersection of the attractor for 
the full dynamics with the invariant subspace, and may be viewed 
as a generalization of on-off intermittency in which the attraction to 
and repulsion from the invariant subspace are controlled by the same 
dynamics. 

Another mechanism for bursting occurs when the strange attrac- 
tor on the invariant manifold attracts typical orbits near the surface 
but is unstable in the sense that there are unstable periodic orbits 
embedded within the chaotic set which are transversely repelling. If 
the trajectory comes near such an unstable periodic orbit, there will 
be a burst away from the invariant surface. Such bursts may occur 
intermittently if noise is present or small changes (called "mismatch") 
are made to the dynamical system that destroy the invariant surface. 
The instability of the first of these orbits is known as a bubbling tran- 
sition (see, e.g., Ashwin et al. |7^, Venkataramani et al. |7^, |80[). 

2.7 Bursts in neural systems 

In neural systems, bursting refers to the switching of an observable 
such as a voltage or chemical concentration between an active state 
characterized by rapid (spike) oscillations and a rest state. Models 
of such bursting typically involve singularly perturbed vector fields 
in which system variables are classified as being "fast" or "slow" 
depending on whether or not they change significantly over the du- 
ration of a single spike. The slow variables may then be thought of 
as slowly varying parameters for the equations describing the fast 



variables (Rinzel ||8l| , Bertram et al. Wang &; Rinzel |84], 
Guckenheimer et al. |2^). As the slow variables evolve, the state of 
the system in the fast variables may change from a stable periodic 
orbit (corresponding to the active state) to a stable fixed point (cor- 
responding to the rest state) and vice versa; such transitions are often 
associated with a region of bistability for the periodic orbit and the 
fixed point, but need not be. Mechanisms by which such transitions 
can occur repeatedly have been classified (Rinzel |]8T| , 82|, Bertram et 
al. [^], Wang & Rinzel |8^]). Behavior of the time interval between 
successive spikes near a transition from the active to the rest state is 
discussed by Guckenheimer et al. |29]; in this paper the presence of a 
subcritical Hopf-homoclinic bifurcation is also identified as a possible 
mechanism for the transition from an active to a rest state. 



3 A new mechanism for bursting 

3.1 Description of the mechanism 

We now describe a bursting mechanism which involves the interaction 
between oscillatory modes in systems with approximate D4 symmetry, 
where D4 is the symmetry group of the square. This mechanism can 
lead to bursts of large dynamic range very close to the instability 
onset (Moehlis & Knobloch |, g, Knobloch & Moehlis g) and is 
expected to be relevant in many different systems with approximate 
D4 symmetry. This symmetry may be present for obvious or subtle 
reasons, as the following discussion demonstrates. 

Consider binary fluid convection in a system of large but finite 
aspect ratio. If the separation ratio is sufficiently negative the sys- 
tem is overstable, i.e., the primary instability is via a Hopf bifurca- 
tion. This is the case for the ^He/^He mixture used by Sullivan & 
Ahlers Q in their experiment carried out in a rectangular container 
D = {x,y,z\ - < X < ^T,-^Ty <y< ^Ty,-^ < z < ^} with 
r = 34, Ty = 6.9. In this experiment Sullivan & Ahlers observed that 
immediately above threshold (e = [Ra — Rac)/Rac = 3 x 10~^) con- 
vective heat transport may take place in a sequence of irregular bursts 
of large dynamic range despite constant heat input (see figure ^) . In 
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Figure 4: Bursts in binary fluid convection with separation ratio 
S = -0.021. (a) e = 3 X lO""^, (b) e = 3.6 x lO'^. The dynamic 
range of the bursts, measured by the range of (A^ — l)/e where N 
is the Nusselt number, decreases with increasing e while the burst 
frequency increases. After Sullivan & Ahlers [4]. Courtesy G. Ahlers. 



this system the presence of sidewalls destroys translation symmetry 
in the x direction which would be present if the system were un- 
bounded, but with identical boundary conditions at the sidewalls the 
system retains a reflection symmetry about x = 0; the primary modes 
are thus either even or odd with respect to this reflection (Dangel- 
mayr & Knobloch [|85|). Numerical simulations of the appropriate 
partial differential equations suggest that the bursts observed in the 
experiments involve the interaction between the first odd and even 
modes of the system (Jacqmin &: Heminger |86]; see also Batiste et 
al. as described in section |3.2D . Thus, to describe the dynamical 
behavior near threshold we suppose that the perturbation from the 
conduction state takes the form 



^{x, y, z, t) = e-2Re {z+{t)f+{x, y, z) + y, z)} + 0(e), (1) 



where e <C 1, f±i-x,y,z) = ±f±{x,y,z). 

Following Landsberg & Knobloch ^ we now derive amplitude 
equations describing the evolution of and Z- using symmetry ar- 
guments. To do this, we first briefly review the topic of bifurcations 
in systems with symmetry (see, e.g., Golubitsky et al. |88| and Craw- 
ford & Knobloch Q). Suppose that 

v = f{v,X), (2) 

where v € i?" and A S R"^ represent dependent variables and system 
parameters, respectively. Let 7 E G describe a linear group action 
on the dependent variables. We say that if 

/(7^;,A)=7/(^^,A) (3) 

for all 7 G G then (^) is equivariant with respect to the group G. 
This is equivalent to the statement that if v(t) is a solution to (^) 
then so is 'yv{t). For example, if a system is equivariant under left- 
right reflections and a right-moving wave exists as a solution, then a 
reflection-related left-moving wave will also exist as a solution. For 
the amplitudes 2;+ and z^ the requirement that a reflected state (ob- 
tained by letting x — > —x in (|^)) also be a state of the system gives 
the requirement that the amplitude equations be equivariant with 
respect to the group action 

Ki : {zj^,Z-) (z+,-z_). (4) 

Moreover, gued by Landsberg & Knobloch |^], the equations 
for the formally infinite system cannot distinguish between the two 
modes, i.e., in this limit the amplitude equations must also be equiv- 
ariant with respect to the group action 

K2 ■■ {z+,z-^) {z_,z+) (5) 



which we call an interchange symmetry. These two operations gen- 
erate together the symmetry group D4. For a container with large 
but finite length, this symmetry will be weakly broken; in particular. 



the even and odd modes typically become unstable at slightly dif- 
ferent Rayleigh numbers and with slightly different frequencies (see 
section ^^). The resulting equations are thus close to those for a 1:1 
resonance, but with a special structure dictated by the proximity to 
D4 symmetry. Finally, we may put the equations for zj^ and Z- into 
normal form by performing a series of near-identity nonlinear trans- 
formations of the dependent variables so as to simplify the equations 
as much as possible (see, e.g., Guckenheimer k, Holmes [^). The 
normal form equations have the additional symmetry (Elphick et al. 



a:[z+,z.)^e'%z+,z^), (TG[0,27r), (6) 

which may be interpreted as a phase shift symmetry. Truncating the 
resulting equations at third order we obtain 

= [\ + i:\\ + i{uj + i:^uj)]z+ + A{\z+\^ + \z-\^)z+ 

+B\z+\^z+ + Cz+zi (7) 

z- = [\- /\\ + i{uj - /\oj)\z. + A{\z+\^ + \z_\^)z- 
+B\z^\^z^+Cz-z\. (8) 

Here Aa; measures the difference in frequency between the two modes 
at onset, and A A measures the difference in their linear growth rates. 
Under appropriate nondegeneracy conditions (which we assume here) 
we may neglect all interchange symmetry-breaking contributions to 
the nonlinear terms. In the following we consider the regime in which 
A, AA, and Aw are all of the same order; in the large aspect ratio 
binary fluid convection context this will occur when these quantities 



are all 0{T~'^) (see section 3^). We will see that when AA and/or 
Au; are nonzero, (0,^) have bursting solutions. Thus, the bursting 
mechanism may be viewed as an interaction between spontaneous 
symmetry breaking (in which the trivial conduction state loses sta- 
bility to a convecting state with less symmetry) and forced symmetry 
breaking (in which the presence of sidewalls make AA and/or Aw 
nonzero, thereby breaking the D4 symmetry). The introduction of 
small symmetry-breaking terms is also responsible for the possibility 



of complex dynamics in other systems that would otherwise behave 
in a regular manner (Dangelmayr & Knobloch pT| , |85[| , Lauterbach 
& Roberts [||, Knobloch (9|], Hirschberg & Knobloch ||). 
To identify the bursts we introduce the change of variables 



z+ = P 2 sm I — I lb — 

^ ^ * 2 4 4 



„i(±</.+V)/2 



and a new time-like variable r defined hy dr/dt = p 
these variables (f7|,|8|) become 



-1 



de 

d(p 
d^ 



In terms of 



-p[2Ar + Br{1 + cos^ 0) + Cr sin^ 6* cos 20] (9) 
-2(A + AA cos 6l)/9^ 

sin 61 [cos 9{-Br + Cr cos 2(p) - Cj sin 20] 
-2AApsin6l (10) 
cos e{Bi - Ci cos 20) - sin 20 + 2Aa; p, (11) 



where A = Ar + iAj, etc. There is also a decoupled equation 
for ilj{t) so that fixed points and periodic solutions of equations (^ 
p^ ) correspond, respectively, to periodic solutions and two-tori in 
equations (0,§). 

In the following we measure the amplitude of the disturbance by 
r = -|- l-Z-P = p ^^] thus p = corresponds to infinite amplitude 
states. Equations (|9[]Tl|) show that the restriction to the invariant 
subspace S = {p = 0} is equivalent to taking AA = Aw = in 
( p!c|JTll ) . The resulting D4-symmetric problem has three generic types 
of fixed points (Swift p5|]): 

• u solutions with cos = 0, cos 20 = 1 

• V solutions with cos = 0, cos 20 = —1 

• w solutions with sin = 0. 

In the binary fiuid context the n, v and w solutions represent mixed 
parity traveling wave states localized near one of the container walls, 
mixed parity chevron (or counterpropagating) states, and pure even 




Figure 5: Examples of (a) u, (b) v, (c) w solutions represented in a 
space-time plot showing the perturbation ^ from the trivial state. 

{9 = 0) or odd {9 = vr) parity chevron states, respectively (Lands- 
berg Sz Knobloch j^). Such states are shown in figure ^ using the 
approximate eigenfunctions 

/±(x) = |e-^^+^^' ± e^^-"j cos ^, (12) 

where 7 = 0.15 -|- 0.025i, L = 80 and < x < ^. Depending on 
the coefficients A, B and C the subspace S may contain additional 
fixed points and/or limit cycles (Swift (9^). In our scenario, a burst 
occurs for A > when a trajectory follows the stable manifold of 
a fixed point (or a limit cycle) Pi E S that is unstable within S. 
The instability within S then kicks the trajectory towards another 
fixed point (or limit cycle) P2 G S. If this point has an unstable p 
eigenvalue the trajectory escapes from S towards a finite amplitude 
{p > 0) state, forming a burst. If AA and/or Aw 7^ this state may 
itself be unstable to perturbations of type Pi and the process then 



repeats. This bursting behavior is thus associated with a codimension 
one heteroclinic cycle between the infinite amplitude solutions Pi and 
P2 (Moehlis & Knobloch [|], Knobloch & Moehhs [|). Examples of 
such cycles are shown in figure ^ Since in such cycles the trajectory 
reaches infinity in finite time the heteroclinic cycle actually describes 
bursts of finite duration (Moehlis & Knobloch IQ]). 

For such a heteroclinic cycle to form it is required that at least 
one of the branches in the D4-symmetric system be subcritical (Pi) 
and one supercritical (P2). Based on the ^He/^He experiments, we 
focus on parameter values for which the u solutions are subcritical 
and the v, w solutions supercritical when AA = Auj = (Moehlis & 
Knobloch Q). When A A and/or Aoj 7^ 0, two types of oscillations in 
{9, (j)) are possible: 

• rotations (see figure 0) 

• librations (see figure |8[) . 

For A > these give rise, under appropriate conditions, to sequences 
of large amplitude bursts arising from repeated excursions towards 
the infinite amplitude {p = 0) u solutions. Irregular bursts are also 
readily generated: figure ^ shows bursts arising from chaotic rota- 
tions. Figure IC provides a partial summary of the different solutions 
of (^-11) and their stability properties; much of the complexity re- 
vealed in these figures is due to the Shil'nikov-like properties of the 
heteroclinic cycle (Moehlis & Knobloch , Knobloch Sz Moehlis ) . 

We now focus on the physical manifestation of the bursts. In 
figure |ll| we show the solutions of figures and |8| in the form of space- 
time plots using the approximate eigenfunctions (0). The bursts in 
figure 11(a) are generated as a result of successive visits to different 
(but symmetry-related) infinite amplitude u solutions, cf. figure in 
figure prKb) the generating trajectory makes repeated visits to the 
same infinite amplitude u solution, cf. figure ^. The former state 
is typical of the blinking state identified in binary fiuid and doubly 
diffusive convection in rectangular containers (Kolodner et al. |9^ , 
Steinberg et al. |^, Predtechensky et al. |Q). It is likely that the 
irregular bursts reported in Sullivan & Ahlers Q are due to such a 
state. The latter is a new state which we call a winking state; winking 



Figure 6: Heteroclinic cycles which exist for AA = 0.03, Ao; = 0.02, 
A = I - 1.5i, B = -2.8 + 5i, and (a,b) C = 1 + i, A = 0.0974, 
(c,d) C = 0.9 + z, A = 0.08461. (a) and (c) are obtained numerically, 
and (b) and (d) sketch the complete heteroclinic networks showing 
all connections. 
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Figure 7: Stable periodic rotations at A = 0.1 for AA = 0.03, Aw 
0.02, A = l- 1.5i, B = -2.8 + 5i, C = 1 + i. 




Figure 8: As for figure but showing stable periodic librations at 
A = 0.1253. 
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Figure 9: Time series and peak-to-peak plot showing bursts from 
chaotic rotations with parameters as for figure |^ but with A = 0.072. 



states may be stable but often coexist with stable chevron-like states 
which are more likely to be observed in experiments in which the 
Rayleigh number is ramped upwards (cf. figure |To|) . For other values 
of AA and Aa; it is also possible to find stable chaotic winking states 
and states which are neither purely blinking nor purely winking (see 
figure 12). 

The bursts described above are the result of oscillations in ampli- 
tude between two modes of opposite parity and "frozen" spatial struc- 
ture. Consequently the above burst mechanism applies in systems in 
which bursts occur very close to threshold. This occurs not only in 
the convection experiments already mentioned but also in the math- 
ematically identical (counterrotating) Taylor-Couette system where 
counterpropagating spiral vortices play the same role as traveling 



waves in convection (Andereck et al. |99|, Pierce & Knobloch [|100(| ). 
In slender systems, such as the convection system described above or 
a long Taylor-Couette apparatus, a large aspect ratio F is required for 
the presence of the approximate D4 symmetry. If the size of the D4 
symmetry-breaking terms AA, Aw is increased too much the bursts 
fade away and are replaced by smaller amplitude, higher frequency 
states (see figure 13). Indeed, if Aw S> AA, averaging eliminates the 
C terms responsible for the bursts (Landsberg & Knobloch Q). From 




Figure 10: Bifurcation diagrams for (a) C = l + z and (b) C = 0.9 + z 
with B, AA, Au; as in figure ^ showing the time-average of r for dif- 
ferent solutions as a function of A. SoUd (dashed) Unes indicate stable 
(unstable) solutions. The branches labeled u, v, w, and qp (quasiperi- 
odic) may be identified in the limit of large |A| with branches in the 
corresponding diagrams when AA = Au; = (insets). All other 
branches correspond to bursting solutions which may be blinking or 
winking states. Circles, squares, and diamonds in the diagram in- 
dicate Hopf, period-doubling, and saddle-node bifurcations, respec- 
tively. 




Figure 11: The perturbation ^ from the trivial state represented 
in a space-time plot showing (a) a periodic blinking state (in which 
successive bursts occur at opposite sides of the container) from the 
trajectory in figure 0, and (b) the periodic winking state (in which 
successive bursts occur at the same side of the container) for the 
trajectory in figure ^. 
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Figure 12: (a) Stable chaotic state with repeated visits to the vicinity 
of the same infinite amplitude state for A = 1 — B = —2.8 + 
5i, C =l + i, AA = 0.03, Auj = -0.02 and A = 0.04, corresponding 
to a chaotic winking state, (b) Stable chaotic state at A = 0.03 
with repeated visits to either the same infinite amplitude state or 
symmetry-related ones. The resulting state is neither purely blinking 
nor purely winking. 



these considerations, we conclude that bursts will not be present if F 
is too small or e too large. However, the mechanism is quite robust 
and even for Auj ^ AA it may still be possible to choose A values so 
that bursts of large dynamic range occur (Moehlis & Knobloch Q). 

It is possible that the burst amplitude can become large enough 
for secondary instabilities not captured by the Ansatz (|l|) to be trig- 
gered. Such instabilities could occur on very different scales and 
result in turbulent rather than just large amplitude bursts. However, 
it should be emphasized that the physical amplitude of the bursts 
is 0(62) and so approaches zero as e j 0, cf. eq. (|^). Thus despite 
their large dynamic range (cf. figure 14), the bursts are fully and 




Figure 13: The perturbation ^ from the trivial state showing stable 
periodic solutions for the parameters of figure |^ except with (a) Am = 
0.1 and (b) Auj = 0.5. From these and figure 11(a) we see that the 
bursts fade away with increasing Ao; and are replaced by smaller 
amplitude, higher frequency states. The amplitude scales are the 
same here as for figure 11(a). 



correctly described by the asymptotic expansion that leads to (0,^. 
In particular, as shown in Moehlis & Knobloch the mechanism is 
robust with respect to the addition of small fifth order terms. How- 
ever, the effects of including D4-symmetry- breaking terms in the cu- 
bic terms in (0j8|) have not been analyzed; these terms can dominate 
the symmetry-breaking terms retained in the linear terms when p is 
small. 
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Figure 14: A very large amplitude burst for the parameters of figure |^ 
except with AA = 0.06, Aw = -0.01. 

3.2 Applicability to binary mixtures 

In view of the motivation for studying systems with approximate 
D4 symmetry described above, it is of interest to examine carefully 
the properties of the linear stability problem for binary fluid con- 
vection in finite containers. In the Boussinesq approximation this 
system is described by the nondimensionalized equations (Clune & 
Knobloch pT| ) 

atu+(u-V)u = -VP + f7i?[0(l + 5) - 5??]z + fiV^u, (13) 
at0+(u-V)e = w^V'^Q, (14) 
9fr/ (u • V)r/ = tS/'^t] ^Sl'^B , (15) 

together with the incompressibility condition 



V • u = 0. 



(16) 



Here u = (u, w) is the velocity field in (x, z) coordinates, P, 9 and C 
denote the departures of the pressure, temperature and concentration 
fields from their conduction profiles, and r] = 6 — C . These equations 
are to be solved in the rectangular domain D = {x,z\ — < x < 

2-L , 2 — ^ — 2 J • 

The system is specified by four dimensionless parameters in ad- 
dition to the aspect ratio T: the separation ratio S, the Prandtl and 
Lewis numbers a, r, and the Rayleigh number R. The boundary con- 
ditions appropriate to the experiments are no-slip everywhere, with 
the temperature fixed at the top and bottom and no sideways heat 
fiux. The final set of boundary conditions is provided by the require- 
ment that there is no mass flux through any of the boundaries. The 
boundary conditions are thus 

u = n • Vr/ = on dD, (17) 

6» = at z = ±1/2, d^9 = 0atx = ±T. (18) 

Here dD denotes the boundary of D. 

Figure |l^ shows the results of solving the linear problem describ- 
ing the stability properties of the conduction state u = 6 = rj = 
for parameter values used by Sullivan & Ahlers |^ in their ^He/^He 
experiment: a = 0.6, r = 0.03, F = 34.0. The flgure shows the neu- 
tral stability curves and corresponding frequencies for the first four 
modes in the range 33.0 < F < 35.0 for S = -0.001 and S = -0.021. 
Observe that when l^l is sufliciently small the flrst two families of 
neutral curves are separated by a gap that is much larger than the 
amplitude of the "braids" within each family (figure p^(a)). This 
is typical of what happens in Rayleigh-Benard convection with non- 
Neumann boundary conditions (Hirschberg & Knobloch [|102| ) and 
makes it easy to justify projecting the fiuid equations onto the first 
two modes that become unstable. However, the situation is not so 
simple. This is because in the case of overstability this behavior does 
not persist for all F or all values of \S\. For larger values of these 
parameters the results take instead the form shown in figures |l^(c,d) 
which show the linear stability results for S = —0.021 and the same 



range of values of F as figures 15(a,b). The modes from the different 




Figure 15: Onset of convection in ^He/'*He mixtures {a = 0.6, r = 
0.03) in large aspect ratio containers, (a) Neutral stability curves and 
(b) corresponding frequencies for the first two even (thick lines) and 
the first two odd (thin lines) modes as a function of the aspect ratio 
r for 5 = -0.001. (c,d) The same but for S = -0.021. Courtesy O. 
Batiste. 



families now cross and the first unstable mode belongs to successively 
higher and higher families when extrapolated to small l^l (Batiste et 
al. [^). Figure |l5|(c) shows the crossing of two even modes involving 
a nonresonant double Hopf bifurcation (figure 15(d)). As discussed 
in detail in Batiste et al. the transition between these two types 
of behavior is mediated by a resonant 1:1 mode crossing at a some- 
what smaller value of l^j. The experimental value of the separation 
ratio from Sullivan &: Ahlers S" = —0.021, therefore corresponds 
to the "crossing" case and the projection of the equations onto two 
modes cannot be rigorously justified except in the neighborhood of 
mode crossing points. 

We denote the growth rates and frequencies of the modes z± 
by X± and lo±. For large aspect ratios, the mode frequencies must 
go like L0± ~ LOoo + ci±T~^ + C2±r~^ + • • •. The fact that the fre- 
quency curves in figure ^(d) are essentially parallel "straight lines" 
implies that 01+ ~ ci_. Therefore, Ai^ = (a;+ — iO-)/2 = 0{T~'^) 
for large F (Batiste et al. |]87|| ). Moreover, as argued in Landsberg &: 
Knobloch IH, the parabolic minimum of the neutral stability curve 
leads to the expectation that AA = ^^(F^^). Thus, in the range 
A = 0(F~^), A, AA, and Alo are all of the same order as F — > 00, 
as required for the applicability of equations (0, ^. Of course, close 
enough to the mode crossing point AA <C Aw, and in this region 
averaging methods can be used to eliminate the {z^z"^, Z-z'^) terms 
from the mode interaction equations (Landsberg & Knobloch |p). 
However, for typical values of AA it appears likely that the system is 
correctly described by equations (^, |^, as hypothesized in Landsberg 
& Knobloch g] and Moehlis & Knobloch |]. 

The first odd and even temperature eigenfunctions for S = —0.021 
are shown in figure 1(: in the form of a space-time diagram, with time 
increasing upward. As in the approximate expression ( [l^ ) the eigen- 
function consists of waves propagating outwards from the center of 
the container. The eigenfunction amplitude has a local minimum at 
the center and increases outwards, peaking near the sidewalls. This 
type of eigenfunction was anticipated by Cross [|103[| and is charac- 
teristic of eigenfunctions in systems with positive group velocity (al- 
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Figure 16: The first even and odd temperature eigenfunctions when 
a = 0.6, r = 0.03, V = 34.0 and S = -0.021 in the form of a 
space-time plot with time increasing upward. Courtesy O. Batiste. 



though, strictly speaking, in a finite system one cannot define a group 
velocity since the allowed wave number is quantized by the sidewalls 
as well as being nonuniform). However, for the present purposes the 
most important observation is that for aspect ratios as large as this, 
the odd and even eigenfunctions are essentially indistinguishable, as 
hypothesized by Landsberg &: Knobloch 

3.3 Other systems with approximate D4 symmetry 

There are a number of other systems of interest where an approximate 
D4 symmetry arises in a natural way and the bursting mechanism 
described in section may be relevant. These include overstable 
convection in small aspect ratio containers with nearly square cross- 
section ( Armbruster []104| , |105| ] ) and more generally any partial differ- 
ential equation on a nearly square domain describing the evolution of 



an oscillatory instability, cf. Ashwin & Mei | 106 |. Other systems in 
which our bursting mechanism might be detected are electrohydro- 
dynamic convection in liquid crystals (Silber et al. |107[| ; T. Peacock, 
private communication), lasers (Feng et al. ]108| ]), spring-supported 
fluid-conveying tubes (Steindl & Troger flOSj ]), and dynamo theories 
of magnetic field generation in the Sun (Knobloch & Landsberg [IIC], 
Knobloch et al. fTTll ). 

Perhaps more interesting is the possibility that large scale spa- 
tial modulation due to distant walls may produce bursting in a fully 
nonlinear state with D4 symmetry undergoing a symmetry-breaking 
Hopf bifurcation. As an example we envisage a steady pattern of fully 
nonlinear two-dimensional rolls. With periodic boundary conditions 
with period four times the basic roll period, the roll pattern has D4 
symmetry since the pattern is preserved under spatial translations 
by 1/4 period and a reflection. If such a pattern undergoes a sec- 
ondary Hopf bifurcation with a spatial Floquet multiplier expi7r/2, 
the Hopf bifurcation breaks D4 symmetry. If the invariance of the 
basic pattern under translations by 1/4 period is only approximate 
(this would be the case if the roll amplitude varied on a slow spatial 
scale) , the D4 symmetry itself would be weakly broken and the new 
mechanism described above could operate. 



Also of interest is the Faraday system in a nearly square con- 
tainer. In this system gravity-capillary waves are excited on the sur- 
face of a viscous fluid by vertical vibration of the container, usually as 
a result of a subharmonic resonance. Simonelli & Gollub |112[ | stud- 
ied the effect of changing the shape of the container from a square to 
a slightly rectangular container, focusing on the (3,2), (2,3) interac- 
tion in this system. These modes are degenerate in a square container 
and only pure and mixed modes were found in this case. In a slightly 
rectangular container the degeneracy between these modes is broken, 
however, and in this case a region of quasiperiodic and chaotic be- 
havior was present near onset. When these oscillations first appear 
they take the form of relaxation oscillations in which the surface of 
the fluid remains flat for a long time before a "large wave grows, 
reaches a maximum, and decays, all in a time short compared with 
the period" . The duration of the spikes is practically independent of 
the forcing amplitude, while the interspike period appears to diverge 
as the forcing amplitude decreases. The spikes themselves possess 
the characteristic asymmetry seen in figures and ^ This behavior 
occurs when the forcing frequency lies below the resonance frequency 
of the square container, i.e., precisely when D4-symmetric problem 
has a subcritical branch. Irregular bursts are also found, depending 
on parameters, but these are distinct from the chaotic states found 
by Nagata [113] far from threshold and present even in a square con- 
tainer. 



4 Discussion 

In this article we have seen that there are many different mechanisms 
responsible for bursting in hydrodynamical systems. The table below 
summarizes the different mechanisms described in terms of properties 
that are most relevant to hydrodynamics. Thus no single mechanism 
can be expected to provide a universal explanation for all observa- 
tions. The bursts found experimentally in Taylor-Couette flow (cf. 
section 2.5) and large aspect ratio binary fluid convection (cf. sec- 
tion 3.1) occur very close to the threshold of a primary instability and 



thus have the greatest potential for a successful dynamical systems 
interpretation of the type emphasized here. We have seen, however, 
that even for fully developed turbulent boundary layers at very large 
Reynolds numbers a dynamical systems approach can be profitable 
(cf. section 2.1). 

Although nearly all of the mechanisms we have described rely 
on the presence of global bifurcations, there are important differ- 
ences among them. For example, the bursts in the wall region of a 
turbulent boundary layer described in section 2.1 are due to a (struc- 
turally stable) heteroclinic cycle connecting fixed points with finite 
amplitude; such a cycle leads to bursts with a limited dynamic range. 
In contrast, in the mechanism of section 3A the dynamic range is un- 
limited. Moreover, the role of the fixed points is different: in the for- 
mer the bursts are associated with the excursions between the fixed 
points while in the latter the bursts are associated with the fixed 
points. Because of the asymptotic stability of the cycle the time be- 
tween successive bursts in the turbulent boundary layer will increase 
without bound unless the stochastic pressure term is included; such 



a stochastic term is not required in the mechanism of section 3.1 



In particular, in this mechanism the duration of the bursts remains 
finite despite the fact that they are associated with a heteroclinic 
connection. This is because of the faster than exponential escape 
to "infinity" that is typical of this mechanism. This is so also for 
the mechanism described in section |2.3| . Both of these mechanisms 
involve global connections to infinity and hence are capable of de- 
scribing bursts of arbitrarily large dynamic range. The models of the 
subcritical transition to turbulence and various types of intermittency 
also produce bursts of finite duration but rely on global reinjection 
which produces bursts of bounded amplitude. 
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Table 1. The different mechanisms discussed in this article. The term 
dynamic range refers to the range of amplitudes during each burst. 
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